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FOREWORD 


The investigation reported herein was conducted at The 
University of Tennessee Space Institute under the sponsorship 
of the National Aeronautics and Space Administration (NASA) 
Marshall Space Flight Center under contract NAS8-33225. The 
NASA Principal Representative for the effort was Mr. David 
Mercier . 

The investigation was conducted by Dr. Virgil K Smith, III 
Principal Investigator using -the computing facilitias of The 
University of Tennessee Space Institute and The University of 
Tennessee at Knoxville. The Principal Investigator was assisted 
in this research by Mr. Muftah Abujelala and Mr. C. Lee Cochran, 
graduate students at UTSI. 

Assistance in providing the basic QNEP engine code and 
answering many questions concerning its modification were prr- 
vided by Mr. Stanley Shapiro and Mr. Mike Caddy of the Naval 
Air Development Center, Warminster, Pennsylvania. 



ORIGINAL PAGE IS 
OF POOR QUALITY 

ABSTRACT 

Advanced horizontal take-off orbital launch vehicles feature 
combined air-breathing and rocket propulsion systems. In order 
to maximize the payload boosted to orbit, an optimization tech- 
nique is required to define the proper engine sequencing over 
the flight trajectory. This research has focused on improving 
the mathematical models of the air-breathing propulsion systems, 
which can be mated with the rocket engine model and ir -'orporated 
in trajectory optimization codes. 

Improved engine simulations provide accurate representation 
of the complex cycles proposed for advanced launch vehicles, 
thereby inci easing the confidence in propellant use and payload 
calculations. The versatile ONEP (Quick Navy Engine Program) 
has been modified to allow treatment of advanced turboaccelerator 
cycles using hydrogen or hydrocarbon fuels and operating in the 
vehicle flowfield. These modifications of the engine code have 
been exercised along with typical installational loss schedules 
to demonstrate the utility of the code for turbofan, augmented 
turbofan and ramjet engine cycles. 

Recommendations have been included for incorporation of 
analytical models of additional, turboaccelerator features in 
the QNEP code. This improved engine code will provide a versatile, 
flexible engine model, both for incorporation in accurate trajectory 
analyses and for assessment of advanced propulsion concepts. 
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INTRODUCTION 


The challenge of transporting large quantities of payload 
into orbit has produced a number of innovative space transporta- 
tion systems [1,2]. Several of the currant concepts for advanced 
orbital launch vehicles feature horizontal take-off, winged 
craft, cowered by combined air-breathing and rocket propulsion 
systems [2,3,4). These combined propulsion systems are used in 

order to exploit the high thrust and specific impulse of air- 

^ » 

breathing systems at low altitudes, while reserving the use of 
the higher thrust, lower specific impulse rocket engines at 
higher altitudes. In order to maximize the payload boosted into 
orbit, optimization techniques must define the proper sequencing. 

Numerous techniques for flight trajectory optimization 
currently exist, including NSEG, A Segmented Mission Analysis 
Program for Low and High Speed Aircraft [5], the Program to 
Optimize Simulated Trajectories (POST), [6], and the Rutowski 
Energy Method [7]. The treatment of the propulsion systems in 
these trajectory calculation methods are limited to scaling of 
tabulated or generalized engine data. Hence, these approaches 
are dependent upon flexible engine codes to accurately treat the 
complex "turboaccelerator" cycles which are common for advanced 
orbital launch vehicles. These composite engines generally 
feature both variable geometry and variable operating modes (8,9, 
10,11,12). The variations in engine geometry are usually found 
in the inlet, ducting, nozzle and compression systems; it is 
likely that the turbine system will also require geometry 
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variability over the wide operating range demanded in turbo- 
accelerator engines. It is also usual for these engines to 
feature several operating modes, including turbomachinery cycles, 
thrust augmentation by duct burning or afterburning, ramjet, re- 
generative cooling, expander cycles, ejectors and ocher exotic 
techniques. It is apparent that turboaccelerator cycles can 
have many possible cycle combinations and much variability. 

. in order to have confidence in trajectory 
optimization involving these propulsion systems, one must use 
engine analytical models with a high degree of fidelity to the 
physical processes and cycle variability being modeled. Only 
then can one expect to adequately define the performance lapse 
rate of the engine, its off-design performance variation, and 
tr.e vehicle/engine installational effects. Since most turbo- 
accelerators utilize hydrogen fuel, it is also important that 
the model incorporate the thermodynamic properties for this fuel. 

A recent analysis of the state-of-the-art in engine model- 
ing techniques and the associated analysis methods [13] revealed 
that two versatile engine modeling codes - NEPCOMP (Navy Engine 
Performance Computer Program) [14] and NNEP (Navy-NASA Engine 
Program) [15] offer the most complete capabilities for treating 
turboaccelerator cycles. Yet, even these current codes do not 
incorporate several of the most important features of the ad- 
vanced air-breathing cycles. As noted earlier, most turbo- 
accelerators utilize hydrogen fuel, along with variable geometry 
inlets and nozzlejj. ,i,n addit_''.on, the engine typically utilizes 
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Thu.. .„ ord« to have a high fidelity sl^latlon of the in- 
= -lled engines, these features fust be treated by the cycle 
^ e . Incorporating these features provides an Improved tool 
or optimising and validating both the propulsion cycle selection 
he propulsion system phasing and. therefore, alloms optimis- 
n P-P^Uant cons^ed and the mass boosted to orbit 


OBJECTIVE 


The Objective of this r'esearch mas to’incorporate in the 
Q^P^engtne code, mhich is a derivative of the versatile KEPCOKP 

analysis ^ — P° detailed, high fldelitv 

analy ts or rnstalled composite air-breathing propulsion svstems ’ 

pec.trcally, the research ctilited the p.HP engine code as the ' 
haselrne simulator and incorporated the following new capabilities 
St cases. (l) extension of the inlet subroutine to cal- 
culate the engine perform.ance when operating i„ the vehicles’ 
uompr^sron and expansion field and applieation of the modified 
program to typical turboaccelerator cycles; (2) identification 

age drag losses and norrle afterbodv drag losses and appli- 
ca ion of these loss schedules to typical turboaccelerator cycles- 
3) ertension of the program to allow simulation of hydrogen' ’ 
and «> applieation of the expanded p.P pi^.^am, ineluding 
2 ydrogen model, the installational losses and the vehiele 
raid effects to typical turboaecelerator building block 
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cycles. The results from each of these areas are discussed in 
the following sections. Additional developments that are re- 
quired to improve the degree of fidelity of the engine models 
and to improve their application have been identified. 
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RESULTS AND DISCUSSION 


This investigation has resulted in improved techniques 
for assessing propulsion system performance of advanced orbital 
launch vehicles. These techniques will greatly assist optimiza- 
tion of the flight trajectories of these craft. 

A. Basic QNEP Model . Initial execution of the basic QNEP 
program on the IBM 360/370 revealed that insufficient data 
fidelity was present in the single precision version to provide 
the required convergence accuracy for design point and off-design 
point engine calculations. This results because the basic word 
length for the IBM unit is 32 bytes, while the CDC 6600 (for 
which the program was originally written) has a 60 byte word 
length. Thus, the CDC 6600, like the Univac 1108, provides 
good data fidelity and sufficient accuracy for the convergence 
logic when configured for single precision execution. 

Following conversion of the basic QNEP code to double pre- 
cision, the engine code was exercised on the IBM 360/370 over a 
wide range of engine design point and off- design point calculations. 
The data cases revealed the same convergence accuracy as the 
single precision versions. 

The basic QNEP program has becT' written for the customary 
engineering units and most calculations have been carried out in 
this system. The results have been converted to SI units, which 
are stated first in this report. The customary engineering units 
are stated afterwards, in parentheses. 
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B. Turboaccelerator Cvcles . There are numerous turbo- 
accelerator cycle configurations which have been proposed for 
application to orbital launch vehicles [3,4,8,9,10,11,12]. In 
most cases, these composite cycles are combinations of the turbo- 
jet, turbofan, augmentor or ramjet cycles which are activated 
sequentially over different portions of the launch and recovery 
mission. The turboaccelerator chosen for analysis with the im- 
proved QNEP model is a turbo ramjet, or, more basically, an 
afterburning mixed-flov? turbofan. A valving system is utilized 
whereby the inlet airflow is directed either totally into the fan 
duct or into both the fan and core flovr passages. For flight 
Mach numbers up to 3.5, the engine operates as a turbofan engine, 
utilizing afterburning for take off and for acceleration up to a 
Mac’' number of 3.5. For subsonic cruise after take o..*f and on the 
return mission, the turbofan operates unaugmented. Above a Mach 
number of 3.5 the engine functions as a subsonic combustion ramjet, 
using the afterburner as the combustor. 

To demonstrate the application of the improved QNEP program 
to these multiple engine operating modes, several "building block" 
engine configurations were exercised using the modified QNEP pro- 
gram. These included design and off-design point performance for 
typical twin-spool turbofans, twin-spool aft fans, augmented tur- 
bofans and ramjets. Schematics of these building block cycles, 
showing the flow station and component identification numbers 
• which were used in the application of QNEP, are displayed in 
Figures 1-3. 
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Figure 1. Engine Schematic - Mixed Flow, Dry Turbofan 



















Figure 2. Engine Schematic, Mixed Flow, Aft Fan, Dry or Augmented. 

























The ease of modifying these building block configurations 
is typified by the augmented aft fan configuration shown in 
Figure 2. Alternate configurations can be constructed by adding 
or subtracting components. For example, from this mixed flow 
turbofan engine, a separate flow turbofan can be configured by 
selecting another nozzle from the ’’spare parts” supply and by 
separating the flow accordingly. 

C. Installational Effects . Exclusion of installational 
effects from the calculation of engine performance can yield un- 
realistic results. The QNEP code was modified to allow inclusion 
of additional installational influences and typical installation 
loss schedules were identified for application to the building 
block cycles. 

1. Vehicle Flowtield Effects . A mathematical model was formu- 
lated and a subroutine constructed to allow the influence of 
vehicle forebody compression and afterbody expansion to be in- 
corporated into QNEP. The model incorporates several operating 
modes since the aircraft/ engine flow field depends upon the 
location of the engine on the vehicle. As shown in Figure 4, 
the major influence of these installational con ;iderations is the 
flow field surrounding the engine inlet and exhaust field. In 
contrast to the operation of the inlet without flow turning (Case 
A) , the turning provided by the fcrebody to the incident super- 
• sonic flow field results in an increase in pressure and density 
downstream of the oblique shock. This flow field surrounds the 
inlet, providing a high pressure and high-density flow field 
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from which the engine draws its mass flow. 

The location of the engine on the vehicle determines the 
extent that the rest of the engine, and in particular, the 
nozzle, is also exposed to the high pressure flow field (Case 
B and C) . The placement of the engine on the mid- to-forward 

of the winged vehicle (Case B) implies that the nozzle 
will be exhausting into the high pressure field; however, an 
engine placement on the aft reaches of the vehicle (Case C) , 
yields a nozzle flow/external flow field inferaction at essentially 
ambient conditions. This ambient field results from flowfield 
expansion over the aft surfaces of the vehicle and engine. 

Therefore, the major objective of the mathematical model of 
the vehicle flowfield is to provide definition of the conditions 
existing downstream of the shock created by the forebody. Exact 
definition of these conditions can be very complex since the 
shock field is dependent upon the specific vehicle geometry. How- 
ever, since typically the major forebody feature preceding the 
engine is the vehicle's wing, it can be assumed that the flow 
turning is essentially two-dimensional (2-D) . Hence, the 2-D 
oblique shock relations are available for treating the flow. 

The usual analysis of 2-D supersonic turning involves the 
use of compressible flow tables and charts (e.g. , NACA 1135). 
However, for the engine model it was deemed desirable to use the 
analytical description of this flowfield, rather than using inter- 
polation of stored data tables. Thus, the equations describing 
an oblique shock field were assembled for coding. In addition, 
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it was desired that only the incident Mach number. and the 
flow turning angle, 6, be specified as inputs, thereby freeing 
the analysis from dependence upon the oblique shock charts for 
determination of the shock angle, 0. However, the equation 
describing the flow field can not be solved explicitly for 0. 


tan 6 


2 cot 0 


' sin^ 0-1] 
M^(y + cos 20) j 


2 


Hence, a numerical solution utilizing the Newton-Raphson method 
was used to solve the relationship for the shock angle based upon 
inputs of Mach number and flow deflection angle [16]. This 
numerical approach has demonstrated excellent stability over the 
Mach number and turning angle range of interest (1.0 < M < 10., 

0 ^ < : < 40 *^) . 

The mathematical model describing the flowfield changes 
across the oblique shock has been incorporated in subroutine 
INLET of the QNEP program. The logic for the modified subroutine 
is shown in Figure 5. If the user desires that the engine per- 
formance calculation incorporate the influence of the vehicle 
flowfield, the flow deflection angle is input in the inlet NAME- 
LIST data by control statement CDAT (15, JCX) . where JCX is the 
component number of the inlet. (If the combination of Mach 
number and flow deflection angle inputs produces a detached 
shock, an error message is generated and the performance calcu- 
lation is made with 6 = 0®.) . If the engine is located such that 
the nozzle exhausts to ambient conditions (Case C) , the control 
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5 . Subroutine INLET Logic Sequen 







st3C6tn6nt CDAT (14.JCX) should be inpu*; es 1. ^ if the nozzle 
exhausts to the forebody shock field (Case B) , CDAT (14.JCX) 
should be input as 0. If the engine is not under the influence 
of a forebody flowfield (Case A). CDAT (15,JCX) input as 0. will 
bypass the forebody effect logic. 

An example of engine performance calculations with the 
forebody flowfield modes and engine location is 
provided in Figures 6 and 7. The variation in thrust and thrust 
specific fuel consumption for Cases A. B and C are displayed over 
a typical altitude and Mach number range for an augmented turbo- 
fan engine. A nominal flow deflection angle of 10° was assumed. 
Both of the performance parameters of these figures reveal the 
significant decrease in performance which results from the engine 
inlet and exhaust nozzle in the shock field (Case B) ; yet. for 
the engine inlet in the shock field and the exhaust nozzle in 

the free stream (Case C) , a significant performance increase 
results . 

Inlet an d Nozzle Afterbody Effects . To enable meaningful 
mission analyses to be conducted on the turboaccelerator engine, 
installed performance for each of the building block engines was 
generated considering the drags associated with a typical inlet 
and afterbody over the flight regime of the engines. Both inlet 
d^3g and nozzle afterbody drag were considered in the 
calculation of the throttle-dependent installation drags associ- 
ated with the engines. Both of these drag results were based 
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upon an isolated nacelle. Any interference effects of nacelle 
placement on throttle-dependent drags were not considered, since 
these effects are a function of a specific vehicle design and 
are beyond the scope of this study. (Upon selection of a given 
vehicle configuration, the integrated inlet and nozzle losses can 
be determined and used in QNEP in the same manner that the iso- 
lated nacelle results are used. Tha data requirements are de- 
fined below.) 

The inlet spillage and nozzle afterbody drags are a strong 
function of the particular vehicle mission-especially the maxi- 
mum Mach number. As evidenced in Figure 8, a more simplified, 
lower design Mach number inlet is less sensitive to inlet-engine 
airflow mismatch during part power subsonic flight [17]. The 
potential losses for a turboaccelerator inlet, (^(Jesign ~ 6-7) 
when operated at off-design conditions can prove disastsrous to 
a vehicle's performance. In addition, higher Mach number vehicles 
require larger exhaust areas to be incorporated into the engine 
nozzles in order to produce maximum thrust. This results in 
relatively more nozzle closure and, hence, boattail drag during 
part power, dry operation. 

Therefore, in collecting parametric data for inlet and 
nozzle performance for the composite propulsion systems which are 
applicable to single-stage- to-orbit vehicles, it is necessary to 
consider inlet and nozzle designs which provide acceptable cycle 
• matching over the entire subsonic- to-hypersonic speed range. 

This eliminates, for example, the application in this study of 
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cal Effect of Inlet Spillage on Subsonic 
Iflc T’-iel Consumption (SFC). 




the vast multitude of data which has teen accumulated for low 
supersonic [18,19], high supersonic. [20,21,22,23,14] and hyper- 
sonic inlet designs [25,26,27,28], However, limited data are 
available from composite propulsion system studies over the 
desired Mach nximber range [9,12,29,30]. 

The operation of a hypersonic inlet at less than its 
maximum mass flow ratio is often required in order to permit 
the air induction system to supply the correct airflow to the 

engine. This may occur because of throttling the engines or 

^ • 

because inlet design requirements conflict for take-off, high 
stibsonic cruise, and supersonic- to-hypersonic acceleration, 
therefore dictating a compromise in the inlet size. The amvount 
of compromise depends entirely upon the particular vehicle, it? 
mission, the engine types and the sophistication of the air in- 
duction system. When estimating the performance of such a 
vehicle operating at less than its supercritical mass flow, the 
inlet spillage drag must be included. 

Spillage drag is defined as additive drag minus cowl suction 
[311. It is a correction that is applied to engine net thrust 
to obtain net propulsive thrust. The cowl suction is largely a 
function of cowl shape for a given mass flow ratio, and acts in 
the thrust direction to cancel some portion of the additive drag. 
Additive drag can be thought of as "subtractive thrust", since it 
must be subtracted from engine net thrust to yield propulsive 
thrust . 

The inlet design selected for this study was chosen from 
three different axisymmetric inlet systems examined by Bencze 
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and Sorensen [29] for app_ication to turboramjet powered hyper- 
sonic cruise vehicles. The inlet data were selected from this 
source following an extensive examination of unclassified 
literature (approximately 200 references). These data were 
chosen because the Mach nijmber range (0 - 6.0) compared closely 
to that proposed for orbital launch vehicles [4]. Also, the 
engine type assumed in the reference closely matches the turbo- 
accelerator cycles considered in the present scudy, - d the 
previous investigation included inlets of various characteristics 
and sophistication. 

The study of the three inlets indent if ied an optimum con- 
figuration based upon the range performance of a hypersonic 
cruise vehicle. The inlet pressure recovery schedule for this 
forward- translating-centerbody design is shown in Figure 9. The 
aiscontinuity in the schedule identifies the starting Mach number 
of the mixed compression inlet system. For Mach numbers greater 
than the starting Mach number, the recovery is the standard 
military specification schedule [32]. 

The schedule of inlet spillage drag coefficient f o * 

the selected inlet as a function of mass-flow ratio and local 
Mach number is shown in Figure 10. The inlet spillage drag is 
non-dimensionalized by the local dynamic head and the inlet 
capture area. The data are based on a combination of experimental 
and theoretical results. The data points indicated by the in- 
dividual symbols and the values at a maos flow ratio of unity 
were provided in the referenced paper. The value at interm.ediate 
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mass flow ratios were approximated based upon trends noted for 
the spillage drag variations of similar inlets [21,22,27,31]. 

The results shov; the usual increase of spillage drag with de- 
creasing mass flow ratio. 

Application of these typical generalized hypersonic inlet 
data to the turboramjet cycles provides a refined estimate of 
propulsion system performance. Selection of a specific inlet/ 
engine configuration for an orbital launch vehicle will provide 
inlet data which can be utilized in QNEP in an identical manner. 

The results of applying these typical hypersonic inlet 
data to the augmented aft fan are shown in Figures 11 and 12. 

The "uninstalled" performance is based upon no spillage drag, 
standard inlet recovery [32], a flow deflection of 10° and ex- 
pansion to ambient conditions. The "installed" performance in- 
corporates the inlet recovery and spillage drag schedules dis- 
cussed above. The impact of incorporating these inlet losses 
is to significantly lower the engine performance; failure to 
incorporate these installation effects can yield a significantly 
optimistic result. 

The installed thrust of a propulsion system is also in- 
fluenced by the difference between the nozzle thrust level and 
the overall afterbody and exhaust system drag. These installa- 
tional effects are in addition to the inlet losses considered 
in the previous section. Ideally, thrust production is a function 
of the nozzle performance parameters, and drag production is a 
function of the aerodynamics of flow over the external body 
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surfaces. In reality, there exists a definite interrelationship 
between thrust and drag of the installed system due to the inter- 
action of the external airs cream with the exhaust stream discharged 
from the nozzle. Thus, the drag of the installed afterbody and 
the exhaust system combination is integrally dictated by both the 
aerodynamic efficiency of the external body geometry and by the 
interference effects caused by the presence of a thrust producing 
exhaust stream. 

primary sources of ^afterbody aerodynamic drag are the 
result of the viscous and inviscid nature of the external air- 
stream over the afterbody surfaces and on any base surfaces [ 31 ]. 
The drag affecting afterbody surfaces can be considered to con- 
sist of skin friction and wave pressure drag components to 
account for the influence of the boundary layer and the potential 
flow field, respectively. The drag affecting, base surfaces is 
most correctly classified as a pressure drag, although the fluid 
mechanisms which produce the pressure forces develop from strong 
viscous interaction of the total flow entering the base region 
from both the afterbody surfaces and the exhaust nozzle. As a 
consequence, base drag is known to depend appreciably upon boundary 
layer characteristics. Then, in order to determine the total 
drag of an afterbody and exhaust system, defined herein in co- 
efficient form as, 

afterbody .^^oattail ^Base 

it is necessary to determine not only the Mach number and static 
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pressure distributions along the entire length of both the 
external body surface and the nozzle wall, but also to have 
knowledge of the manner rn which the boundary layer develops 
along each Oi. these surfaces as well. Each drag component is 
referenced to the local dynamic head and the maximum cross- 
sectional area of the afterbody. 

Because of the mission and cycle variations of composite 
propulsion systems, the nozzle systems of these devices will 
require large variations in pressure ratio and, hence, will re- 
quire capabilities for large area ratios [9.30,33,34,35]. This 
trend is illustrated in Figure 13. Yet, large nozzle exit areas 
will result in large drag increments [30,36,37,38]. This requires 
a compromise between m.aximum propulsion efficiency and engine 
installation drag. In order to provide this compromise, it is 
likely that the turboaccelerator nozzle that will be utilized 
in isolated nacelles will be a variable geometry design of either 
the plug or expansion-deflection type [10,11,33.39,40,41]. With 
the incorporation of variable geometry, these types of nozzles 
provide both control of the nozzle throat area (through a trans- 
lating and/or collapsing centerbody) and the throat- to- exit area 
ratio, while retaining a fixed shroud and nozzle exit diameter. 
These features provide the large nozzle area variations required 
of composite propulsion systems over the subsonic to hypersonic 
speed range. 

However, the data on these nozzle types which is available 
in the open literature is very limited. The results are generally 
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confined to performance measurements in the static mode and in 
the transonic- to- low supersonic flight regimes [42 - 51). The 
dearth of high velocity results in the literature likely stems 
from the lack of specific applications in this flight regime 
and the potential of these nozzle tjrpes for engine infrared 
radiation suppression. This lack of information is unfortunate 
since, as noted earlier, the nozzle off-design performance and 
drag is very sensitive to the vehicle's maximum Mach 
number, and, hence, the nozzle pressure ratio and geometry. This 
restricts the usefulness of nozzle results cited above because 
these data relate to the lower velocity range with low nozzle 
pressure ratios and limited variable geometry. Alternately, 
analytical methods for treating the range of typical turboacceler- 
ator nozzle conditions are restricted by the lack of emoirical 
correlations [52,53,54]. 

However, noting these limitations, the nozzle data docu- 
mented in Reference 49 were selected for application to the tur- 
boaccelerator cycles. The selected data is for a plug nozzle 
incorporating both internal and external expansion, having a 
design total- to-static pressure ratio of 25. It is likely that 
the turboaccelerator cycles will have a maximum pressure ratio 
on order of magnitude larger than this; however, variable geometry 
features in the turboaccelerator nozzle will allow the unit to 
configure for optimum performance for a large range of pressure 
ratios. In addition, the limited data on plug nozzles indicate 
that the major influence of external effects occur in the transonic 
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range. This is an area of nozzle performance which will re- 
quire detailed test and analysis for j election of an optimized 
turboac-elerator configuration. 

The selected nozzle data have been extended over the 
pressure ratio range of 4 to 24 using trends of similar data 
from References 43, 49 and 50. These data extensions are indi- 
cated by the dashed portions of the curves in Figure 14. 

The afterbody drag characteristics for the selected nozzle 
are shown in the upper portion of Figure 14, depicting the varia- 
tion of afterbody drag coefficient as a fv -tion of nozzle pressure 
ratio and the external flowfield Mach number. The data exhibit 
a typical rise in afterbody drag at transonic speeds and the 
typical decrease in drag as the design pressure ratio is achieved. 
In some cases, this decrease in drag can result in "negative drag", 
or thrust, resulting from significant forces acting on the base 
areas of the plug nozzle. 

The thrust coefficient data for the selected nozzle as a 
function of nozzle total pressure ratio are shown in the lower 
portion of Figure 14. The thrust coefficient is the ratio of 
the actual nozzle thrust to the ideal thrust of the nozzle flow. 

The ideal thrust equals the actual mass flow rate times the ideal 
velocity (i.e., the velocity reached by the stream upon isentropic 
expansion from the total pressure to the ambient pressure.) The 
data reveal the typical plug nozzle behavior - a cusp in the 
thrust coefficient values which is associated with the region of 
nozzle pressure ratios where the relative significance of internal 
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expansion increases. Qualitatively, the plug nozzle displays 
performance which is superior to the convergent- divergent nozzle 
at all pressure ratios and provides the altitude compensation 
that is desired in a multi-mode nozzle. The data have been 
assumed to be invariant with nozzle throat-to-exit ratio, though 
QKEP can accept data in this format without program modification. 

These internal and external nozzle performance data have 
been incorporated as data tables in the specifying data of the 

typical engine cycles. The results of applying these nozzle in- 

^ * 

stallation effects are typified by the augmented aft fan results 
shown in Figures 15 and 16. The performance is again compared 
on the basis of uninstalled and installed performance. Unin- 
stalled data are based upon no spillage drag, standard inlet 
recovery [32] , a flow deflection of 10*^ and expansion to ambient 
conditions assuming a constant thr^-^st coefficient of 0.99. The 
installed performance incorporates the thrust coefficient and 
afterbody drag schedules discussed above. The impact of incor- 
porating these nozzle losses is to lower the engine performance 
in the region of low nozzle pressure ratio. At higher pressu’-e 
ratios (Ref. Figure 14), the negative drag results in an increase 
in engine performance. Again, it should be noted that these re- 
sults reveal the importance of incorporating individual component 
losses in the cycle calculations. 

The results of applying both the inlet and nozzle loss 
schedules to the augmented aft fan cycle are shown in Figures 17 
and 18. The uninstalled performance is based upon no spillage 
drag, standard inlet recovery 1321, a flov: deflection of 10° and 
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expansion to ambient conditions assuming a constant thrust co- 
efficient of 0.99. The installed performance incorporates the 
inlet recovery, spillage drag, thrust coefficient and afterbody 
drag schedules. The optimistic performance provided by uninstalled 
assumptions is obvious. 

The effect of these combined losses at part-power conditions 
is also significant. Figure 19 displays the installed and unin- 
stalled performance for throttled operation of the unaugmented 

turbofan at M = 0.8 and an altitude of 10973 meters (36,000 feet). 

^ • 

Again, large losses are demonstrated. 

These results have used typical inlet and nozzle loss 
schedules chosen from the open literature to demonstrate the 
importance of incorporating their influence in calculations of 
engine performance. Selection of a specific inlet/engine/nozzle 
afterbody configuration for an orbital launch vehicle will provide 
inlet and nozzle data which can be utilized in QNEP in an identi- 
cal manner. 

D. Hydrogen Combustion Model . Because of the high energy 
content of hydrogen, its lar ;2 heat capacity and its utility for 
both airbreathing and rocket cycles, most turboaccelerator cycle 
concepts utilize hydrogen as the fuel (8,9,10,11,12]. For this 
reason, the QNEP thermodynamic model has been modirfied to allow 
the user to make engine calculations with hydrogen fuel. The 
model has also been expanded to allow calculations with hydro- 
carbon fuels with a broad range of user-specified carbon-to- 
hydrogen ratios. The evaluation of various potential hydrogen 
combustion moaeit comprised a major portion of this effort. The 
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factors which were considered in this evaluation are detailed 
in the paragraphs below. 

The operating environment of turboaccelerators implies 
both increased cycle temperatures and operating Mach numbers 
over those currently encountered by gas turbine engines . At 
these operating conditions, losses of 20 - 307, of the net thrust 
can occur due to the energy absorbed in dissociation. Hence, 
it becomes extremely important chat the effects of chemical 
dissociation are incorporated in the cycle calculation. As 
the degree of dissociation is dependent upon the temperature 
and pressure of the combustion gases, the determination of the 
equilibrium composition, the mean molecular weight and the re- 
sulting thermodvnamic properties must include consideration of 
this dependence. In addition, it is likely that the fuel v.^ill 
be utilized to cool the engine internal structure (turbine blades , 
duct walls, nozzle walls, etc.) and the external surfaces of the 
vehicle. During this use, the enthalpy of the fuel is increased, 
and it is desirable to incorporate provision in the thermodynamic 
model to treat this increased energy content. 

The standard QKEP program provides thermodynamic calcula- 
tions via FUNCTIOK THERM. The method is based on separate cal- 
culation of the thermodynamic properties of air and the products 
of combustion of a stoichiometric fuel-air mixture having the 
hydrogen- to-carbon ratio (0.16) that is typical of jet fuels; it 
is then asstaned that the properties of the combustion products 
for any fuel-air ratio that is less than the stoichiometric 
value may be obtained by a linear interpolation between the tnco 
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extreme values. In the absence of dissociation, the results of 
this approach are highly accurate; hcv^ever, application of this 
approach to cycles in which dissociation occurs will give only 
approximate results. Since the high operating temperatures of 
advanced turboaccelerator cycles will likely result in signifi- 
cant dissociation, large errors in engine calculated performance 
can result. For this reason, FUNCTION THERM is restricted to 
hydrocarbon fuel used for cycle temperatures below 4500*^F, and 
processes where significant dissociation is. absent. Hence, its 
calculated thermodynamic properties are functions only of tem- 
perature and fuel-air ratio. This factor led to the conclusion 
that, for the range of temperatures and pressures encountered 
by turboaccelerators, it would not be possible to merely sub- 
stitute the hydrogen combustion product properties in the current 
FUNCTION THERJI logic. Instead, it was necessary to modify the 
complete program logic for determining the thermodynamic prop- 
erties. 

A number of methods for computing thermodynamic properties 
of exhaust gases are available for computer use. They fall into 
two general categories : one contains tabulated combustion prop- 

erties, while the other is a computer program which calculates 
the properties. Examples of the former data format include 
Keenan and Kaye [55], Banes, et al. [56], Powell, et al. [57], 
and Browne and Warlick [58]. Examples of the latter include 
Pinkel and Turner [59], Osgerb.y and Rhodes [60], Gordon and 
McBride [61], Pelton [62], and Mascitti [63]. The use of the 



Cabular form of pre-calculated daca using some form of curve 
fitting was rejected because of the large data storage require- 
ments; it is anticipated that approximately 4000 constants would 
be required in o-der to provide temperature, pressure and fuel- 
air ratio dependence for hydrogen fuel, sacrificing both computer 
storage and accuracy. 

Combining the QNEP code with a highly generalized code like 
Gordon and McBride [61] could greatly hinder its usefulness be- 
cause of the significant increase both in computer storage and 
operating requirements and in execution time. The basic double 
precision QNEP code requires approximately 256 kilobytes of 
storage for execution; incorporation of the NASA SP-273 code 
in the QNEP code would require approximately 500 kilobytes of 
3-orage ^or execution. In addition, the thermodynamic subroutine 
of QNEP is a function which is interrogated by every component 
subroutine in the program. Hence, it is used more often and can 
be more time consuming than any other subroutine in the engine 
model. Thus, its computational efficiency is a significant factor 
in program utility. Calculation of the thermodynamic properties 
at single values of pressure, temperature and fuel-air ratio using 
the thermodynamic model of Reference 61 requires approximately 
one second. In considering the numerous iterations of thermo- 
dynamic properties required in cycle calculations, this computa- 
tional requirement is considered too excessive. Additionally, 

• the Gordon and McBride code contains many capabilities which are 
not required for engine cycle calculations. Elimination of these 
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excess cape.bilicies could prove very complex and time consuming. 
In contrast, the incorporation of an engine-specialized computer 
code for calculating the thermodynamic properties provides im- 
proved computational efficiency and convenience, and, in addi- 

» retains the tlexibility of orderly updating as compositions 
for fuel or air change, or as thermodynamic property data are 
improved. 

These factors have led to the selection of the Mascitti 
model [63] for incorpor ation ^in the v^NEP engine program. The 
code is a simplified combustion gas model which includes the 
®^fficts of dissociation and which allows a wide range of fuel-air 
ratios and carbon- to-hydrogen ratios. The model allows treatment 
of hydrocarbon-air combustion, hydrogen-air combustion and 
dissociating air. The model is applicable over the pressure 
range of 0.001 - 100 atmospheres and up to 7000°F. The gas model 
is simplified by neglecting the formation of species containing 
atomic nitrogen, thereby enabling a considerable simplification 
of the composition equations and allowing a solution for the 
chemical composition to be obtained with a single-level iteration 
of these equations. This assumption was made since calculations 
for this temperature and pressure range indicate that the forma- 
tion of nitrogen species, such as N, NH, and NO, occurs in 

negligible amounts, and, therefore, has a very small effect on 
the thermodynamic properties of combustion gas mixtures [8,9]. 

The impact of this assumption was examined both -in terms of the 
thermodynamic properties of stoichiometric kerosene-air and 
hydrogen-air combustion products and in terms of the performance 
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of an idealized subsonic c'^mbustion ramjet. Good agreement was 
demonstrated between the simplfied model and more comprehensive 
treatments in the range of temperatures applicable to hypersonic 
engine cycles [63]. 

Preliminary analysis of the model revealed that several 
errors in constituent properties were present in the program 
listing. These errors were eliminated using data from the lates 
JANAF tables [64]. In addition, routines for calculation of 
specific heat properties vere added to the program, and entry 
routines to allow the program to make various cycle calculations 
were constructed and evaluated. Since the calculated thermo- 
dynamic properties are functions of temperature, pressure, fuel- 
air ratio and carbon- to-hydroger. ratio, each of the com.ponent 
sutroutines required modirication to allow interfacing with 
the thermodynamic lunction subroutine. Tliis new thermodynamic 
subroutine, which was designated FUNCTION THERMO, demonstrated 
an average execution time of 0.07 seconds for calculation of the 
thermo dynam.ic properties at single values of pressure, tempera- 
ture and fuel-air ratio. 

The results of these modifications to the QNEP code were 
evaluated by comparing typical engine cycle calculations using 
both calculator-based design computations and results from, the 
standard QNEP engine code. Because FUNCTION THERMO has the 
capability for treating both hydrocarbon and hydrogen fuels , 
it was possible to execute typical engine cycles with hydro- 
carbon fuel (carbon- to-hydrogen ratio of 0.16) for the modified 
QNEP program and compare the resulting data with the standard 


}r h 


QNEP program (containing FUNCTION THERM) . For cycle conditions 
for which no dissociation occurs, comparison of the outputs from 
the two programs provide an excellent check of the validity of 
the program modifications and the accuracy of the simplified 
model incorporated in THERMO. 

Typical results of cycle calculation comparisions tor an 
augmented aft fan engine are shown in Table 1. These data 
tabulate the percentage of error encountered in thrust and thrust 
specific fuel consumption for calculations using the QNEP code 
incorporating the two thermodynamic routines. The values of 
error are given by: 


Error = 


THERM Value-THERMO Value 
THERM Value 


(For simplicity, the thrust is tabulated in pounds, while the 
thrust specific fuel consumption is cited in pounds per hour 
per pound.) Excellent agreement is demonstrated between the 
values calcuxa^ed from the two models, with a maximum discrepancy 
in thrust of appro cimately Z%. The maxim.um discrepancy in thrust' 
specific fuel consumption is less than 1%. 

By artificially restricting the dissociation calculation 
in FUNCTION THERMO, it v?as also possible to obtain data for 
comparable high altitude and Mach number conditions with the ram- 
jet cycle. These data w’hich are displayed in Table 2 also re- 
veal excellent agreement for the performance parameters. 

Comparable calculations were conducted for the front fan, 
unaugmented turbofan cycle for both maximum power and throttled 
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Table 1. Comparison of Thermodynamic Models for 

Augmented Aft Fan Performance Calculation 

ALTITUDE = 0.0 m. (0.0 ft.) 



THERM. 

THERMO . 

7o ERROR 1 

MACH 

FN 

TSFC 

— 

FN 

TSFC 

FN 

TSFC 

o 

o 

102901 

1.625 

102256. 

1.636 

.6268 

-.6769 

.05 

101617 

1.65 

100964. 

1.661 

.6426 

-.6667 

.15 

99861 

1.696 

99182. 

1.708 

.6799 

-.7075 

.25 

99763 

1.735 

98986. 

1.747 

• .7788 

-.6916 

.35 

101912 

1.759 

100681. 

1.776 

1.2079 

-.9665 

.45 

105144 

1.779 

103991. 

1.794 

1.0966 

-.8432 

.55 

109105 

1.796 

108394. 

1.808 

.6516 

-.6682 


ALTITUDE = 3048 m. (10,000 ft.) 



THERM. 

THERMO . 

% ERROR 

MACH 

FN 

TSFC 

FN 

TSFC 

FN 

TSFC 

.4 

76124 

1.761 

75514 

i.772 

.801 

-.624 

.6 

81664 

1.8 

80776 

1.813 

1.087 

-.666 

.8 

94897 

1.788 

93344 

1.805 

1.636 

- .9507 

1.0 

104067 

1.818 

102980 

1.833 

1.04 

-.769 

1.2 

113181 

1.856 

117756 

1.868 

.817 

- .646 

1.4 

121679 

1.9 

120541 

1.913 

.935 

-.684 

1.5 

127899 

1.916 

126571 

1.928 

1.03 

-.626 
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Table 1 (Continued) 
ALTITUDE = 9144 n. (30,000 ft.) 



THERM. 

THERMO. 

% ERROR 

MACH 

rN 

TSFC 

FN 

TSFC 

FN 

TSFC 

.6 

41455 

1.766 

41026 

1.778 

1.034 

-.679 

.9 

51859 

1.754 

51157 

1.766 

1.353 

-.684 

1. 

56747 

1.745 

55819 

1.757 

1.635 

-.687 

1.2 

68074 

1.729 

66516 

1.745 

'2.288 

-.925 

1.4 

79408 

1.729 

77885 

1.742 

1 . 917 

- .751 

1.6 

87865 

1.754 

86362 

1.766 

1.71 

-.684 

1.8 

95589 

1.789 

i 

93940 

1.8 

1.725 

-.614 


ALTITUDE = 15,240 m. (50,000 ft.) 



THERM. 

1 

THERMO . 



\ ERROR 1 

MACH 

FN 

TSFC 

FN 

TSFC 

FN 

TSFC j 

.8 

19243 

1.754U 

19119 

— 

1.763 

.644 

1 

- .398 

1.1 

25313 

1.724 

25020 

1.734 

1.157 

-.58 

1.6 

38638 

1.708 

37998 

1.719 

1.656 

- . 644 

2.2 

52463 

1.785 

51552 

1.796 

1.736 

-.616 

2.6 

61845 

1.856 

61436 

1.861 

.661 

-.269 

3. 

80276 

1.893 

78101 

1.90.*^ 

2.709 

-.633 

3.2 

94423 

1.898 

91602 

1.912 

2.987 

1 

-.737 


























Table 2. Comparison of Thermodynamic Models for 
Ramjet Performance Calculation 

ALTITUDE = 15240 m. (50,000 ft.) 



THERM . 

THEPJ^O. 

% ERROR 

M.A.CH 

FN 

lOgil 

FN 


FN 

TSFC 

3.0 

6109^ 

2 .275 

60741 

2.279 

.5778 

- .1758 

3.5 

110135 

2.209 

109652 

2.22 

.43855 

-.498 

4. 

182485 

2.186 

181699 

2.196 

.4307 

- .4575 

5. 

3S5113 

2 . 232 

3S3804 

2.24 

• .3399 

-.3584 

6. 

550176 

2.414 

547475 

2.425 

.491 

-.4557 

7. 

442300 

3.055 

438807 

3.082 

.7897 

- .884 


ALTITUDE = 21,366 m. (70,000 ft.) 



THERM. 

THEPJMO . 

% ERROR 

MACH 


TSFC 

FN 

TSFC 

FN 

TSFC 

3.0 

23295 

2.278 

23140 

2.282 

.6654 

-.1756 

3.5 

41973 

2.212 

41754 

2.223 

.522 

- .497 

4. 

69514 

2.189 

69149 

2.199 

. 525 

-.457 

5. 

146427 

2.236 

145802 

2.245 

.427 

- .4025 

6. 

208174 

2.422 

206959 

2.434 

.58365 

- .4955 

7. 

— J 

163970 

i: 

3.094 

162472 

3. 122 

.9136 

' .905 

























Table 2 (Continued) 


ALTITUDE = 27,432 m. (90,000 ft.) 



THERM. 

THER.MO . 

7, ERROR 

M.\CH 

FN 

TSFC 

FN 

TSFC 

FN 

TSFC 

3. 

8934 

2.292 

8854 

2.298 

.8955 

-.262 

3.5 

16060 

2.226 

15942 

2.238 

. 7347 

- . 539 

4. 

26500 

2.204 

26342 

2.215 

.746 

-.499 

5 . 

55418 

2.256 

55052 

2.266 

* .6604 

.443 

6. 

76969 

2.460 

76325 

2.473 

.837 

-.5285 

7. 

54520 

3.301 

53795 

3.339 

1.33 

-1.1512 


ALTITUDE == 30,480 m. (100,030 ft.) 
THERM. THERMO. 7, ERROR 






























power conditions. These data comparisons also revealed a 
maxitnxmi discrepancy in thrust of 3% for maximum power conditions. 
For throttling of the engine, this discrepancy typically in- 
creased, with a maximum discrepancy of 57. in thrust observed 
near the lOTo thrust throttle setting. Comparable errors in 
thrust specific fuel consumption were also observed. 

The excellent agreemenc between the THERNO and THERM models 
applied in QNEP for a hydrocarbon fuel gave high confidence for the 
validity of the program modifications and the accuracy of the simpli- 
fied model incorporated in THERMO. In addition, the execution of 
the QNEP code for hydrogen fuel was evaluated by comparing the 
results from the code with calculator-based cycle computations. 
Again, excellent agreement was demonstrated. In addition, Mascitti 
(63] shov;ed excellent agreement with the tabulated results of 
Bro\im and Warlick (581 in the usual operating ranges of turbo - 
accelerator cycles. 

The previous experience with the code [12,63] and the 
author’s evaluation led to the conclusion that the model was in- • 
deed performing well and that application to typical building 
block cycles should proceed. 

E. Modified QNEP Code Application . Following verification 
of the THERMO subroutine and the individual thermodynamic modifi- 
cations to the QNEP subroutines , the installational modifications 
and the thermodynamic modifications were combined. After check- 
out of the combination, the modified QNEP code, which is documented 
in Appendix B, was applied to typical turboaccelerr ♦•or building 
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block cycles. These cycles, which are depicted schematically in 
Figures z and 3, are tj'pical turbofau, augmented turbofan and 
ramjet cycles Each of the selected cycles are described by 
their QNEP data inputs in Appendix C. The modifications to these 
data inputs, which are required of a user in order to treat in- 
stallational effects end to allow treatment of various fuels are 
described in Appendix A. Appendix C also contains sample output 

for each cycle type illustrating the use of the modified QNEP 
program. 

Typical results from these applicati.ons of the modified 
QNEP code are depicted in Figures 20-26. The losses due to in- 
corporating installational losses are evident in each figure, in- 
^be e*.^ect on part— pov.’er nerformance shor-.C; in Figure 
These results again confirm the importance of incorporating 
reaj. world installational losses in estimates of engine per- 
formance. 

Figures 25 and 26 also display the importance of incorporat- 
ing real gas losses in the engine model. The data depicted in 
symbols illustrate tne additional losses i.n performance due to 
dissociation of the combustion gas products and resizing of the 
exhaust nozzle to handle the equilibrium flow. The cases labeled 
installed and uninstalled are based upon artificially precluding 
dissociation in the thermodynamic model and sizing the nozzle to 
handle the frozen flow. (This mechanism is controlled by varying 
in FUNCTION THERMO the temperature level considered to be the 
boundary for the onset of dissociation effects). Isentropic ex- 
pansion to ambient conditions is assumed in both cases. 
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The QNEP code offers great versatility and flexibility. 
The modified QNEP code extends its canabilities and utility for 
treating turboaccelerator cycles. 
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CONCLUSIONS AND RECOMMENDATIONS 

The versatile ONEP code has been modified to give greater 
utility for treating advanced turboaccelerator concepts therebv 
pro\ading an enhanced tool for conducting trajectory and payload 
optimization calculations. The inclusion of hydrogen fuel pro- 
vides a major step toward exact modeling of turboaccelerator 
cycles. The capability to incorporate the effect of vehicle 
.-lowfield on the engine perrormance has been demonstrated to be 
significant for accurate cyc'-e calculations! Finally, the in- 
clusion Oi. typical inlet and nozzle losses has been somewhat 
hindered by the lack of specifically applicable data in the 
open literature; however, selection and application of repre- 
sentative losses has demonstrated the importance of identifvins 
uheir contribution to the overall propulsion system analysis. 

It is important that other turboaccelerator features be 
modeled and incorporated in QNEP in order to increase confidence 
in the engine calculations . Chief among these currently neglected 
features are regeneratively cooled nozzles, pre-compressive cooling 
and turbo-expander energy extraction. In view of the importance 
of advanced orbital launch techniques to the further utilization 
Oj. the orbital environment, it is recommended that timelv develop- 
ment of the additional analysis capabilities be pursued. 


60 


' 4 . 


REFERENCES 


1. Love, E. S. , "Advanced Technology and the Space Shuttle," 

AIAA Paper No. 73-31, January 1973. 

2. Shelton, B. W. , "Advanced Earth- to-Orbit Transportation 
Vehicles and Their Propulsion Systems," AIAA 78-1102, pre- 
sented at the AIAA/SAE 14th Joint Propulsion Conference, 

Las Vegas, Nevada, July 25-27, 1978. 

3. Jackson, L. R. , J.A. Martin, and W. J. Small, "A Fully 
Reusable, Horizontal Takeoff Space Transport Concept with 
Two Small Turbojet Boosters,” NASA TbI 74087, October 1977. 

4. Hanley, G, M. , "Future Space Transportation Systems," 
presented at the Canaveral Council of Technical Societies, 
Fifteenth Space Congress, Cocoa Beach, Florida, Apr:’ 1 26-28, 
1978. 

5. Hague, D. S. and H. L. Rozendaal, "NSEG-A Segmented Mission 
Analysis Program for Low and High-Speed Aircraft," Volume 
1-Theoretical Development, NASA CR-2807; Volume 2-Program 
Users Manual, NASA CR-2808; August 1977. 

6. Brauer, G. L., D. E. Cornick, and R. Stevenson, "Capabilities 
and Anolications of the Program to Optimize Simulated Tra- 
jectories (POST), "NASA CR-2770, February 1977. 

7. Rutowski, E. S. , "Energy Approach to the General Aircraft 
Performance Problem," J. of the Aeronautical Sciences, March 
1954, pp. 187-195. 

8. Bendot , J. G. , P. N. Brown, and T. G. Piercy, "Composite 
Engines for Application to a Single-State- to-Orbit Vehicle," 
NASA CR02613,* December 1975. 

9. "A Study of Composite Propulsion Systems for Advanced Launch 
Vehicle .Applications," Vols. 3 and 6, Marquardt Report No. 
25194, September 1966. 

10. "Estimat d Performance of a Mach 8.0 Hydrogen Fueled Turbo- 
fanramjet (STFRJ-230A) " Pratt & Whitney Aircraft Report, No. 
TDM-1905, 1-25-65. 

11. "Supersonic - Hypersonic Propulsion System Study Data, Tandem 
Turboramjet (GeV/JZI - Study B) ," General Electric Company, 
SS-63-5, April 15, 1963. 

12. Waters, M. H. , "Turberamjet Propulsion System for All-Body 
Hypersonic Aircraft," NASA TN D-5993, January 1971. 


61 



% 


. ^ t 

^ I 


13. Smith. V. K. , "Air Breathing Engine/Rocket Trajectory Opti- 
mi?.ation,” in "Faculty Fellow Research Reports, 1978 NASA/ASEE 
Sunmer Faculty Fellow’ship Program," University of Alabama BER 
Report No. 224-94, September 1978. 

14. Shapiro, S. R. and M. J. Caddy. "NEPCOMP - The Navy Engine 
Performance Program," ASME Paper 74-GT-83, March-April 1974. 

1.5. Fishbach, L. 11. and M. J. Caddy. "NNEP - The Navy NASA 
Engine Program.," NASA TM X- 7 185 7, December 1975. 

16. Zucrow, M. and J. Huffman, Gas Dvnam.ics , John Wiley 6t Sons, 
1976. 

17. Johnson. J. E. . "Variable Cycle Engines - The Next Step in 
Propulsion Evolutiou?," AIAA Paper /6-758, 1976. 

• 

18. Petersen, M. W. and G. C.* Tamplin, "Experimental Review of 
Transonic Spillage Drag of Rectangular Inlets," USAF Report 
APL-TR-66-30, May 1966. 

19. Hawkins, J. E. , F. P. Kirkland and R. L. Turner. "Inlet 
Spillage Drag Tests and Numerical Flow-Field Analysis at^ 
Subsonic and Transonic Speeds of a 1/8- Scale, Two-Dimensional, 
External-Compression, Variable-Geometry, Supersonic Inlet 
Configuration," N.AS.A. CR-268C, .\pril 1976. 

20. Wasserbauer, J. F. and W. H. G<=>rs tenmaier , "Inlet-Engine 
Matching for SCAR Including Application of A Bicone Variable 
Geometry Inlet," AIAA Paper 78-961, July 1978. 

21. Woolett, R. R. , E. T. Meleason, and D. A. Choby, "Transonic 
Off-Design Drag and Performance of Three Mixed Com.pression 
Axisj’xnmetric Inlets," NASA TM X-3215, June 1975. 

22. Woolett. R. R. . E. T. Meleason, and D. A. Choby, "Transonic 
Off-Design Drag and Performance of an Axisymmetric Inlet With 
40 Percent Internal Contraction on Design," NASA TM X-3042, 
August 1974. 

23. Smeltzer, D. B. and N. E. Sorenson. "Tests of a Mixed Com- 
pression Axisvmmetric Inlet With Large Transonic Flow at 
Mach Number O'. 6 to 2.65," NASA TN D-6971. December 1972. 

24. Smeltzer, D. B. and N. E. Sorenson, "Analytic and Experimental 
Performance of Two Isentropic Axisvmmetric Inlets at Mach 
Number 0.8 to 2.65," NASA TN D-7302, June 1973. 

25. Pearson, L. W. and L. F. -Jillv, "Hypersonic Research Engine 
Project. Phase 2A, Inlet Program," Final Technical Data 
Report, 1 March 1967 - 28 February 1969, NASA CR 66797, 

March 1969. 


62 


26 . 


•r.. " • ' '• 

, r 

• • *■ » ^ ' V 4. • * 


27. 

28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38 


Trexler, C. A., "Inlet Performance of the Integrated 
Langley Scramjet Module," AIAA Paper 75-1212. ^ 

^^ssell, E. A. Macklev. and A. L 

ESe^A4-otL™H" "he NASA Hypersonic Research 

L?emLr 1974 Integration Model." NASA TM X-3038, 

Bontrager, "Wind Tunnel Tests of A 
8,"'lESclR-69-9,'Fe^^^^ A. 5.6. and 

Me?it^of ■Th;^^? "A : udy of the Relative 

•Mlssiln?" MlI%aperlS!6l7' Cruise 

VphSff’ ^^^'^ere. and R. Ceresuela. "Propulsion of 

Vehicules Hypersoniques . "^ AGARD LS-42, Vol. 2, July 1972. 

Rockweil^CorDoraM^r''?i/^v°'^^^^^^ Manual." North American 
Julv IQftR u?’ <^o^^hus Division, Report NR68H-434, 

July 1968, (Available as AD 837501). 

Tu?bof^i ^’^glnes. Aircraft. Turbojet and 

u bofan. Military Specification Mil-E-5009D. November 1967. 

Hanley G. A. and R, Burgeron, ".\n Overview o'^^ the Satellit-A 

AlAA7sAE^14?h’jo-^rp^^®"r:’’ '5-975, presented 

NevkSi, m- 27 19?r Conference, Las Vegas, 

> Hypersonic Ramjet Development " in Conbus t-’ on 

mian: 

^’^hbage, W. A. Vahl, and U J Small 
^Studies of Airframe-Propulsion-Svstem Integrccion ^or Mach 6 
^ruise Vehicles." NASA TN D-4128, October jf67. ^ ^ ^ 

Probfomc'^A ?•• "hypersonic Aircraft and Their Aerodvnamic 
roblems. in Progress in Aeronautical Science^; Vo1 fi 
Pergamon PressTTybb. pp. 271-353'; ~ 

Duper. G. L. "Comparison of Hypersonic Ramjet Engines wirv, 
Subsonic and Supersonic Combustion." in Combustion^nd Pro^^ 

f^nr^A. Milan r9‘6U, Pergamon P«ss. 


r^^TT~pP- 84-119 
Johnston. P. J. . j. h. Cubbage. and J. P. Weidner "Studi-s 

AIAA^apor No Atrckafr," 


« 


63 


39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48. 

49. 

50. 


Bendot J G P. N. Brown and T. G. Piercy, "Cotnposite 
Ingines for Application to a Single-Stage-to-Orbit Vehicle. 

NASA CR-2613, December 1975. 

Calmcn. J.. C. Menioux, 

du Turbostatoreacteur , in AGAKD ur ^ i, i^oo. hf 

Chansler, R. C. "Hypersonic "’'naIaI 

Phase 2A. Nozzle Program. Terminal Summary Report, 

CR-101532, December 1968. 

Staid. P. S.. "Wind Tunnel Performance Tests of Coannular 
Plug Nozzles," NASA CR 2990. April 1978. 

Wasko R A "Performance of Annular Plug and Expansion- 
Deflection nAzziL Including External Flow Effects at 
Transonic Mach Numbers," NASA TN D-4462^, Apr^.1 196 . 

Hoad V L "Fliaht Investigation of Installation Effects 

Nozzle Si ^h a Serils of Boattailed Primary Snrouds 
llstahld o^L UndCrw^ing Nacelle." NASA TM X-2983, February 
1974. 

Harrin--on D. E. . S. M. Nosek. and D. M. Straight." Cold-Flow 
Perfor^-ance of Several Variations ^ '' 

Nozzle for Supersonic-Cruise Aircraft, i^ASa i.. X ..ixU, 

October 197A. 

Berrier B L and C. E. Mercer. "Off-Design Performance 
of Two isentroCic Plug Nozzles Designed for a Pressure Ratio 
of 16.5," NASA TN D-3852, March 1967. 

... T 3 A n F HAr'rin&ton *’Perforrnance of a Collapsible 

?!:rAozzlt Hfving- E?thef i.«?Po;itlon Cylindrical or.,Var-,ab.e 
Angle Floating Shrouds at Mach Numbers from 0 .o 2.0, a ^ 

TM X-1657, October 1968. 

Mercer C B. and L. B. Salters, "Performance of a Plug Nozzle 
Ravine’ a Concave Central Base With and Without Terminal 
rllSngt a? T?aLon?c Speeds." NASA TH D-180A. May 1963. 

Valerino A S , R. F. Zappa, and K. L. Abdalla, "Effects of 
External ’stream on the Performance of 

Nozzles at Mach Numbers of 2.0. 1.8. and 1.5. NASA ME.iU 
20-17-59E, March 1959. 

Hearth D P. and G, C. Gorton. "Investigation of IHrust and 
Drarcharacteristics of a Plug-Type Exhaust Nozzle, NASA 
RM E53L16. February 1954. 


6A 


. ' . T 

„ > *X . A X 1 


51. 


52. 


53. 

54. 


55. 

56. 


57. 


58. 


59. 

60. 


61. 


62. 


63 . 


Stoddard J. A. P., "Jet Effects on Boattail Pressure Drag 
at Supersonic Speeds," in AGARD-CP-91-71 , ’Inlets and Nozzles 
for Aerospace Engines," December 1971, pp. 23/1-23/11. 

Glasgow, E. R. , J. S. Divita, P. C. Everling and J. A. Laughrey , 
"Analytical and Experimental Evaluation of Performance Pre- 
diction Methods Applicable to Exhaust Nozzl^s.^ Paper 
No. 71-719, June 1971. ■ 

Migdal, D. , "Supersonic Annular Nozzles," AIAA Paper No. 

71-43, January 1971. 


Gaede, A. E. and J. J. Lopez, "Selection of Nozzle Contours 
for a Research Scramjet Engine," AIAA Paper No. 67-453, 

July 1967. 

Keenan, Joseph H. , and J. Kaye, Gas Table s, John Wiley and 
Sons, Inc., 1948. 

Banes, B. , R. W. McIntyre, and J. A. Sims, "Properties of 
Air and Combustion Products with Kerosine and Hydrogen 
Fuels," Bristol Siddeley Engines Limited, 1967. 


Powell, H. N., and S. N. Suciu, E. , Properties of Combustion 
Gases/System; C^H 2 ^-Air. McGraw-Hill Company . Inc., 19o5. 


Bro'.cn , W. G. and D. L. 1-,'arlick, 
Gases/System; H.,-Air," General 
1962. 


"Prooerties of Combustion 
Electric Report R62FPD-366, 


Pinkel, B. and R. L. Turner, "Thermodynamic Data for the 
Computation of the Performance of Exhaust-Gas Turbines , 

ARR No. 4B25, October 1945. 

Osgerby, I. T. and R. P. Rhodes, "An Efficient Numerical 
Method for the Calculation of Chemical Equilibrium in the 
H/C/O/N/A System," AEDC-TR-71- 256 , April 1972. 

Gordon, S. and B. J. McBride, "Computer Program for Calcula- 
tion of Complex Chemical Equilibrium. Compositions, Rocket 
Performance, Incident and P^eflected Shocks, and Chapman- 
Jouget Detonation," NASA SP 273, 1971. 

Pelton, J. M. , "A Computer Model for the Calculation of . 
Thermodynamic Properties of Working Fluids of A Gas Turbine 
Engine," AEDC-TR-76-15, March 1976. 

Mascitti, V. R. , "A Simplified Equilibrium Hydrocarbon-Air 
Combustion Gas Model for Use in Air-Breathing Engine Cycle 
Computer Programs," NASA'TN D-4747, 1968. 


k. 


65 


Stull, D. R. and H. Prophet, JANAF Thermochemical Tables , 
Second Edition, U. S. Government Printing Office, June 1971. 

Caddy. M. J. and S. R. Shapiro. "NEPCOMP-The^ N.avyp»^^ 5 ie 
Performance Computer Program, Version I , ” ao 
740A5-30, April 1975. P P'-A'R 



^CEDIHG PAGE 


3UM( NOT famed 


APPENDICES 




APPENDIX A 

MODIFIED QNEP ENGINE CODE 
USERS' GUIDE 


I 


USERS' GUIDE 


The increased capabilities of the modified QNEP engine 
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APPENDIX C 

SA:!?L£ input and output 

Cl. UNAUGMENTED AFT FAN TURBCFAN 
C2. AUGMENTED AFT FAN ENGINE 
C3. RAMJET ENGINE 
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ENGINE PERFORMANCE CALCULATION 
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C2. AUGMENTED AFT FAN TURBOFAN 

SAMPLE INPUT AND OUTPUT 
FOR THE 

ENGINE PERFORMANCE CALCULATION 
WITH INSTALLATION EFFECTS AND H 2 FUEL 
(DISSOCIATION ARTIFICIALLY PRECLUDED) 
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COAT(l,21 1*0.,0.,0.,1.,0.,H,, 
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FlCSiTsl. «iNO 
AFT FAN tNGlNt, DCSIC4 POINT 
50 FLO* SCHEOUbC 

I » 0. 

T 1 


• # 1 • » 
piooO 


io.. 


# 1 • ^ 0 0 


THCT 21 


WRA 21 


COT 

75 

Z 1 
NACH 17 


HCOR 1 


0. 

• 5 

.75 

1.157 

1.06 

1.054 

.911 


.4 

.10 

1.55 

1.06 

1.046 

.796 


RAW RECOVCRT 

0 . 

0 . 


1.75 

4.5 

2000 . 


.25 

1.9S 

5.0 


.57 

.675 

.66 

• 665 

#•5 

.90 

• 95 

.975 

1.54 

1.735 

1.93 

2.12 

1.06 

1.06 

1.06 

1 .06 

1.041 

1.033 

1.022 

1 .014 

.70* 

.644 

.596 * 

• 556 

.5 

.75 

1.0 

1.25 

2.0 

5.2 

2.5 

3.0 

3.5 


«69. 

1.0 

2.31 

1.06 

1.0 

.520 


1.5 

4.0 


• • 


RR 

1 

l»0 

RR 

1 

1.0 

RR 

1 

1*0 

RR 

X 

1.0 

RR 

X 

1.0 

RR 

X 

.01 

RR 

X 

.04 

RR 

X 

••• 

RR 

X 

.12 

RR 

1 

.R2S 

RR 

1 

.17 

RR 

X 

.•1 

RR 

X 

.74 

RR 

X 

.67 

RR 

X 


RR 

X 

Isi 

RR 

COT 

X 

.40 


100 

1 

MACH 


iT ' * *• 

•• I. M ^ * 


SPILLAGE DRAG 

0. 


RPR I 

t'OS i 

RPR 7 
CDS 7 
RFR 7 
CDS 7 
RFR 7 
CDS 7 
•FH 6 
COS 6 
MPR 4 
CDS 4 
MFR 3 
CDS 3 
MFR 4 
CDS 4 
MFR 1 
CDS 1 
EOT 
120 

Z 1 


• 3S 
1. 

• 4 

• 5 

• 6 

.7 

.8 

.0275 

0. 

• 022 

• 014 

• 0095 

.006 

.0035 

. 4 

.5 

• 6 

.7 

• 8 


.045 

• 0265 

.017 

.0105 

.006 

.0002 

. 4 

• 5 

• 6 


• 8 


.0S4S 

• 033 

.022 

• 014 

.008 

.0003 

. 4 

• 5 

• 6 

.7 

.8 

• 9 

.072 

• 045 

.0285 

• 016 

• 0105 

.0005 

• s 

• 6 

.7 


.9 

1 , 

.057 

.035 

• 022 

• 012 

• 0005 

0. 

• 7 

• 8 

.9 

1. 



• 028 

• 013 

.0005 

0— - 



.8 

• 9 

1. 




• 009 

.0002 

0. 




.8 

• 9 

.973 

1» 



.017 

« 

• 0003 

0. 

0. 




0 . 

NOZZLE THRUST COEFFICIENT 

0, 


T 

1 

0. 



PR 

13 

3. 

4. 

5. 

cv 

13 

14. 

16. 

18. 

• 94 

• 91 

.89 



.97 

.979 

• 983 

COT 

.. 




150 

AFTERBOOT 

DRAG 


z 

1 

0. 


MACH 

3 

1.5 

1.8 

2.0 

NPR 

11 

4« 

6. 

8. 



18. 

20. 

22. 

COAB 

11 

• 131 

.125 

.107 



• 02 

• 015 

• Oil 

NPR 

10 

6. 

8. 

10. 



20. 

22. 

24. 

CDA8 

10 

.079 

.07 

.056 



• 0 

-.003 

-.005 

NPR 

11 

5. 

6. 

8. 



18. 

20. 

22. 

COAB 

11 

• 061 

.059 

• 052 



-.009 

-.019 

-.017 


SO. 

24« 

.04 

-.010 


czz 








• • • • 


• • • • 

I* K* K , , 

COX 





200 

« VS. 

CURBKCTCD 

rLO«»fAN u:>fcD 

FOR bPRsl*2 

Z X 

0. 




Y 1 

1. 




H 2 

0. 

.5 

1.1 


WDA 1 

0. 

.5 

1.1 

* 

LOT 





20S 

B VS. 

erncitNCY 

•fan ds^d fur 

BPksi.2 

1 i 

0. 




Y 1 

1. 




R id 

.4 

.5 

.6 

.7 .1 .9 .929 


.952 

.0$ 

.V9 

1. 1.1 

err 12 

.93: 

.842 

.851 

.ab» .198 ,89« .198 


• 8565 

• 8548 

• 853 

.89 .884 

eoT 





300 

B VS. 

CO»«ECTPD 

FLO- COMPRESSOR 

z t 

0. 




Y i 

1. 




k 3 

0. 

.5 

1.011 


kkk 3 

0. 

.5 

1.011 



eoT 


305 

Z 

Y 

1 

1 

R VS. 
0. 

1. 

EFFECIENCY-COMPRESSOH 





R 

9 

.416 

l.O 

•522 .633 

1.011 

.76 

• 889 

.944 

• 968 

trr 

9 

• 666 
• 81 

.74 ,774 

• 804 

.794 

• 811 

.813 

• 811 


EOT 


CNO OF TABLE RATA 

AO ALTsO.,t,**sl,rtPal .TDiCLsO. .Xi'AsO.O 
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1. irn,i‘i^ = t.,r, iTt;.j*>) = i.,Cw«JU.n)si., 

Cu*Ill.leJ*l..Ci'4ru,l-)al..Cr.ilU,iO) = l..CiJAU'. A»».C 

Ar I PA.4 UiSIG'. P"I4T. CONlP'JLS Ot* 

40 Mpsl,TneLsO.,TiT2 3l$O..ALTsO..COA1 ( 4 , li )s32»0. , 

TlIL£*l,lPKll«rsO, 

CDATO.l 3)s.9b,C0AT(A. 13 )s 49900.,CDAT(10,U)*1000.,CDAT(U ,l3)=0.0,»tB=l2, 
AHAsO.,.Os, .1 ,.15..2..2S,.3,.35..4. .4S..5..55 AfcttO 
AiJUxeuTED AFT FAN 

40 ALTs20O0.,M*!=12,X*<As.2..3,.4..5,.6,.7,.B,.9,l.,l.i,1.2,1.3 4tN0 
40 ALTS4000. 4fe.N0 
40 ALIS60OQ. 4 EnD 
40 AblsfeOOO. nKKiO 

40 ALT=iOOOU.,Xy4:.4..S..4..7..|i..9,l.,l.l,1.2.1.3, 1.4,1.5 4CMD 
40 ALT=20000.,X-4 As. 6,.R,.9,1. ,1.1 ,1.2, 1.3, 1.4,1. b, 1.6,1. 7,1. • 4tND 
40 ALTs30000.,A^A=.h,.8..9,l..l.l,1.2,1.4,1.6,1.7,1.8,2.0,2.2 4felN0 
40 ALrs40000.,XMAs.B,.9.1.,l.l,1.2,1.4,1.4,2.2,2.4,2.b,2.6,2.8, 4END 
40 ALl^bOOOO. ,XfiA3.3,.9,1.1.1.4,1.6,1.6,2.2,2.4,2.6,2.»,3.,3.2 4END 
40 ALT360000. 4EN0 

40 AblsTOOOO. ,XMA:.e,l.,1.2,1.6,2.0,2.2.2.4.2.ii.3.0.3.2,3.4.3.S 
4EM0 

40 ALl-SSOOO. 4END 
40 bNOKUUsi 4bN0 
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